


SUBTASK REPORT Po 


UGE A ASE 


00067 3260 


PROBLEMS AND POSSIBILITIES OF 
HIGH ALTITUDE WORLDWIDE MAPPING 


16 JANUARY 1967 


<>” SIMULATED ERROR ANALYSIS: 2 o-qe 
<0 °. UTILIZING MATHEMATICAL = 7S ty 
SYSTEM GEOMETRY MODELS — 





Declassified and Released by theN RO 





in Accordance with E. O. 12958 © 


TAR erence nerr sO, BREE 
5a NOV 2 61997 


—FOP-SEGRE RUF 


SUBTASK REPORT 
Copy Ne 


PROBLEMS AND POSSIBILITIES OF 
HIGH ALTITUDE WORLDWIDE MAPPING 


16 JANUARY 1967 


SIMULATED ERROR ANALYSIS 
UTILIZING MATHEMATICAL 
SYSTEM GEOMETRY MODELS 


Itek 


LEXINGTON 73, MASSACHUSETTS 


ae a 
FANDCE“vrie— 


FOP-SEPRE EE CONTROTSYS#EM-~ONLY 


ITEK CORPORATION 








1. INTRODUCTION 


On 10 December 1965, Subtask Report for Task 1, Subtask C, Development of Mathematical 
System Geometry Model as completed and submitted. The document 
explained the philosophy e mathematical exercises required in the generation of error 


analysis programs for determining the mapping potential of calibrated panoramic (PG) and 
metric imagery. ; 


The error analysis programs, based on the theoretical panoramic calibration, were 
completed by June 1966, and simulated tests were run at that time. Prior to June, however, the 
error analysis programs were expanded to provide a mathematical tool for performing control 
intensification in the operational mapping task. 


Subsequently, due to technical calibration complications, the philosophy of panoramic 
calibration was revised, thus requiring reprogramming of the system geometry models to make 
them operationally applicable. However, the error analysis, based on the theoretical calibra- 
tion concept, still provides a useful tool for pan-metric systems analysis. 


The purpose of this report is to present an example of the use of the system geometry model 
computer program, as an error analysis tool. In the particular example shown here, the cap- 
ability of calibrated panoramic materials to produce accurate relative point locations is estimated 
and compared with a given set of map accuracy specifications. 


The system selected is a J-1 PG unit with internal geometry calibrated to +4 arc-seconds 
(1-sigma level) flown at an operational altitude of 85 nautical miles. All other mission and 
system parameters are in accordance with real operational data. 
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2. SYSTEMS GEOMETRY MODEL 


The pan-metric adjustment of the systems geometry model is a classic weighted block 
adjustment modified to functionally constrain the camera stations to a simple Keplerian orbit 
and to allow inclusion of convergent panoramic photography. Interior orientation elements for 
both the frame and panoramic cameras are included as adjustment parameters (as are the 
exterior orientation angles of both cameras) with those for the panoramic cameras expressed 
as cubic functions of time. Time is included as an independent measured variable for the 
panoramic model; i.e., each image of interest on a panoramic photograph will not only have 
a measured x and y coordinate, but also a measured time coordinate. 


From the standpoint of describing the physical conditions existing during a panoramic 
camera exposure, i.e., the position and orientation of the camera and ground coordinates 
changing as functions of time, the adjustment is unique. The application of this adjustment 
technique is at present limited to the number of panoramic photographs which lie within one 
metric model due to computer storage limitations and the fact that Keplerian elements will 
adequately approximate the true position of the camera as a function of time only for short 
periods. 


The required inputs to this adjustment are photocoordinates and the variances of points 
of interest (photocoordinates of images on panoramic photographs are defined as including a 
time coordinate as well as x and y coordinates). Estimates of all parameter values are also 
required. As desired, variances and covariances of all parameters can be included as input. 


In general, the output from this adjustment is the covariance matrix of the adjusted 
parameters; however, due to the data processing technique utilized, the contribution to ground 
positioning error from the photographic measurements and the parameter can be—and is— 
separated. All internal computations take place in a siderial geocentric system, but for 
interpretative reasons ground coordinates with their attendant variance-covariance matrices 
are transformed into three other systems: (1) geocentric x, y, and z; (2) geographic latitude, 
longitude, and elevation; and (3) cross-track, in-track, and elevation. 


This program has been mor rjbed in Task 1-C, Development of Mathematical 
Systems Geometry Models, December 1965. 


The following tables list the various input and output conditions for the pan-metric adjust- 
ment program. 


Table 1 presents the assumed standard error to be expected in each of the variables entering 
the model. The program utilizes these errors to provide constant weighting factors in the least 
squares fit. , 


Two covariance matrices [Tables 2(a) and 2(b)| are presented in terms of Keplerian elements 
as well as cross-track, in-track, and elevation standard deviations. The first of these matrices 
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is the input covariance matrix obtained from the tracking data or orbit determination program. 
The matrices listed in Table 2(a) were obtained from the Aerospace Trace Program. 


Since the pan-metric program is an adjustment program, the possibility of improvement 
of orbital elements exists, based upon the photogrammetric condition equations. The second 
covariance matrix {Table 2(b)| is the covariance matrix of the adjusted orbital elements. The 
improvement can be noted by comparison of the diagonal elements. 


Table 3 presents the standard deviation for the ground point locations to be expected from 
the reduction. Table 4 presents the errors in distance and elevation from the centroid of the 
ground point distribution tothe point as functions of the distance from the centroid. 





Table 1 — Input Covariance Data 


J-1 PG 
Orbit 85 nm 
Photo points frame 
Oxi = Ovi = 25 x 107° mm* 
Photo points pan 
ot; = 0%) = 169 x 107° mm* 
Frame camera orientation (local system) 
On = Oo = 10 arc-sec, 
Pan camera orientation (local system) 


Unconstrained 


Parameters held to constant o for all cases 


on = 9x 10-® sec? 


Times oR = 25 x 107!" sec? 
oT, T, = 0.899873 x 107° 
o?, = 2.7976 x 10°" rad’ 
Ground 42 = 5.5573 x 107!” rad? 
points a 
oj, = 92903.04 m* 
(A o of 106.5 meters corresponds to a 90 percent 
Xp o=5x 107? mm, 
Frame Yp a9=5x 107? mm, 
F > =15 x 107? mm, 
PXp os =1x 1078 mm 
PYy a =5%x 10-> mm. 
Fan PF o = 1.5 x 107? mm, 
Camo = 2 = 107? mm, 





(o = 5 x 1072 mm) 


13 x 107? mm) 


(o = 


Gy) = 35 arc-sec 


(o = 3 ms) 
(o = 50 psec) 


(Covariance to express relative time) 


(o = 106.5 m) 

(a = 106.5 m) 

(o = 1000 ft) 
C.E.P. of 750 feet) 
o” = 25 x 1078 

o = 25 «x 1078 

o” = 225 x 1078 

go” = 25 x 1078 

ot = 225 x 1078 

o = 4x 1078 
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Table 3 — Output Standard Deviations 


Ground Someters meters ‘meters 
Point No. Cross-Track iIn-Track Elevation 


12 Reid cece 4%,7 
i¢ 1. Z1.? a%,a 
<A tm La? 47.% 
2? wi,’ 47.7 ai ek 
ae w1,? er. 44.7 
$e Sobek diet 4i+7 
Sn ra 276% au. 
37 %1.7 24.4 40.5 
as &1,1 O40 4764 
64 14. edet 4s 
6? Sra, 14.7 45.4 
79 44.1 des 470s 
74 ad? 13.3 4204 
7S we L303 atek 
7% tat que AN, 
32 Pear zie” ac0,5 
Qh ier 274% a°.5 
9e tease S$ z4iet at,4 
94 obeys 2404 4105 
LiF Lie? one ak. 
120 19.9 13.7 43.0 
124 Ve 17.2 42.9 
126 Wile 4 17.1 41.7 
130 Va. Lien 41.1 
13? apa riaes 17. 45.9 
134 aie ae Zte7 4766 
13¢ 5 a ook ar eR 
142 fo hg} 2s. 4°.5 
144 5 oe ete a',& 
174 aera alee Que 
174 oh yt 17,6 47,48 
272 ua pera 47.4 
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194 aaa ete 6747 
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22g ae Par 6%,7 
727 VW gs 4g 42747 
234 TAs Jie aes 
23h rd aie 13,- 4444 ‘ 
24% Lperd gles 4164 
24% ee” ete 4°,9 
744 M1. cis Gi, 
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25? Sg! cle at. 
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Table 4 — Distance and Elevation Standard Deviations as a 
Function of Distance from Centroid 


Ground Ometers Miles Ometers 
Point No. Distance Distance Elevation 
13 tw." 7%.7 nD] 
3A $265 75. BN .D 
BA Bi fees Sieh 4784 
«48 Soe Pe 5447 
224 Be ©. 130° Seda 
144 phet Biies 6908 
179 125% 3n.7 5: a 
228 33.° 7467 beet 
25? 33° 77,2 50,3 
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